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Abstract —We propose a secure transmission scheme for a 
relay wiretap channel, where a source communicates with a 
destination via a decode-and-forward relay in the presence of 
spatially random-distrihuted eavesdroppers. We assume that the 
source is equipped with multiple antennas, whereas the relay, 
the destination, and the eavesdroppers are equipped with a 
single antenna each. In the proposed scheme, in addition to 
information signals, the source transmits artificial noise signals 
in order to confuse the eavesdroppers. With the target of 
maximizing the secrecy throughput of the relay wiretap channel, 
we derive a closed-form expression for the transmission outage 
prohahility and an easy-to-compute expression for the secrecy 
outage prohahility. Using these expressions, we determine the 
optimal power allocation factor and wiretap code rates that 
guarantee the maximum secrecy throughput, while satisfying a 
secrecy outage prohahility constraint. Furthermore, we examine 
the impact of source antenna number on the secrecy throughput, 
showing that adding extra transmit antennas at the source brings 
about a significant increase in the secrecy throughput. 

I. Introduction 

Wireless communications are inherently insecure, due to 
the broadcast nature of the medium, which makes security 
a pivotal design issue in the implementation and operation 
of current and future wireless networks. Compared to the 
traditional key-based cryptographic techniques that are applied 
to upper layers, physical layer security can enhance the secrecy 
of wireless communications without using secret keys and 
complex encryption/decryption algorithms, and thus has been 
recognized as an alternative for cryptographic techniques. The 
key idea of physical layer security is to exploit the randomness 
of wireless channels to offer secure data transmissions [1, 
2]. In early studies, e.g., [3], the principle of physical layer 
security was established in a single-input single-output wire¬ 
tap channel. Subsequently, physical layer security in multi¬ 
input multi-output (MIMO) communication systems has been 
intensively investigated [4-15], due to the benehts of MIMO 
techniques such as high data rate and high link reliability. 

Most recently, physical layer security in large-scale wireless 
networks, such as mobile ad hoc and sensor networks, has 
been receiving considerable attention [16-19]. A key property 
of large-scale wireless networks is that the node locations 
in the network were spatially randomly distributed. As such, 
stochastic geometry and random geometric graphs are used 
for modeling the locations of spatially random-distributed 
nodes. In [16], the throughput of large-scale decentralized 
wireless networks with physical layer security constraints was 


investigated. Considering the path loss as the sole factor 
influencing the received signal-to-noise ratio (SNR) at the 
receiver, [17] examined the secrecy rate in cellular networks. 
In [18] the secrecy rate achieved by linear precoding was ana¬ 
lyzed in the broadcast channel with spatially random external 
eavesdroppers, and in [19] the secrecy rate achieved by linear 
precoding in cellular networks was analyzed. 

We note that [16-19] only considered point-to-point trans¬ 
missions. This leaves physical layer security with cooperative 
relays in large-scale networks as an open problem. Since the 
relay channel efficiently improves the coverage and reliability 
in wireless networks [20,21], it is of practical signihcance 
to investigate the secrecy performance of such channels. 
Particularly important in this context would be extensions of 
previous work on relay wiretap channels which focused on 
scenarios where the location of an eavesdropper(s) is hxed 
and known at the source (e.g., [22-27]). 

In this work we propose a secure transmission scheme for 
a relay wiretap channel, where a source transmits to a desti¬ 
nation via a decode-and-forward (DE) relay in the presence 
of multiple spatially random-distributed eavesdroppers. The 
source is equipped with multiple antennas, while the relay, 
the destination, and the eavesdropper are equipped with a 
single antenna each. In addition to information signals, we 
assume that the source transmit artihcial noise signals in 
order to confuse the eavesdroppers. Aiming at maximizing 
the secrecy throughput, while satisfying a secrecy outage 
probability constraint, we determine both the optimal power 
allocation factor between AN signals and information signals 
and the optimal wiretap code rates. 

Compared to current studies on physical layer security in 
relay networks [22-27], our contributions are threefold. Eirst, 
we consider a more practical scenario where the locations 
of the eavesdroppers are spatially randomly distributed and 
not known at the source. Second, we propose a new secure 
transmission scheme that maximizes the secrecy throughput 
in such a scenario. Third, we derive explicit expressions for 
the transmission outage probability and the secrecy outage 
probability, which enable us to analytically characterize the 
secrecy throughput of the relay wiretap channel. We note that 
these expressions are independent of realizations of the main 
channel and the eavesdropper channel. 

The rest of the paper is organized as follows. Section 
II describes the system model and presents the proposed 


o 


o 


o 



o O 

o 



Fig. 1; Illustration of a relay wiretap channel with spatially 
random eavesdroppers. 


scheme. In Section III, we analyze and optimize the secrecy 
performance achieved by the proposed scheme. Numerical 
results and related discussions are presented in Section IV. 
Finally, Section V draws conclusions. 

II. System Model AND Proposed Scheme 

We consider a relay wiretap channel, as depicted in Fig. 1, 
where a source (S) communicates with a destination (D) with 
the aid of a relay (R) in the presence of multiple spatially 
random eavesdroppers. We assume that the source is equipped 
with N antennas, while the destination, the relay, and the 
eavesdroppers are equipped with a single antenna each. We 
also assume that there is no direct link between the source and 
the destination. We denote dsr and drd as the distance between 
the source and the relay and the distance between the relay 
and the destination, respectively, and denote rj as the path loss 
exponent. As in [16-19], the eavesdroppers are modeled as 
a homogeneous Poisson Point Process (PPP) $ with density 
A. This model is practical and representative for decentralized 
networks where each node is randomly distributed [28]. The 
source, the relay, and the destination do not belong to 

A. Transmission Scheme 

We now detail the proposed transmission scheme. The 
proposed scheme utilizes two time slots. In the first time slot, 
the source transmits information signals together with AN 
signals to the relay. In this time slot, the transmitted signals 
from the source are overheard by the eavesdroppers. The role 
of the AN signals is to confuse the eavesdroppers. In the 
second time slot, we assume that the source is silent. While it 
is clear that this will lead to a sub-optimal solution, such an 
assumption does allow for analytical tractability in a special 
case. Removing this assumption means more power can be 
added to the signal in the hrst time slot. Our other work [29], 
in which we consider multiple antennae at all nodes, quantihes 
the advantage of adding source noise in the second time slot. 
We also assume the relay transmits the received signals to 
the destination using the DF protocol. In this time slot, the 
broadcast signals are also overheard by the eavesdroppers. 
We assume that all the channels are subject to independent 
and identically distributed Rayleigh fading. We also assume a 
quasi-static block fading environment in which all the channel 


coefficients remain the same within one time slot. We denote 
hsr as the 1 x TV channel vector from the source to the relay 
and denote h^d as the channel coefficient from the relay to the 
destination. 

We now express the transmitted and received signals in two 
time slots separately. In the first time slot, the transmitted 
signal at the source is given by 

xs = Wt, (1) 

where W denotes the NxN beamforming matrix at the source 
and t denotes the combination of the information signal and 
the AN signal. To perform a such transmission, we hrst design 
W as 


W = [ws Wan] , (2) 


where ws is used to transmit the information signal and Wan 
is used to transmit the AN signal. The aim of W is to degrade 
the quality of the received signals at the eavesdroppers. By 
transmitting AN signals into the null space of through 
W, the source ensures that the quality of the received signals 
at the relay is free from AN interference. In designing W, 
we choose ws as the eigenvector corresponding to the largest 
eigenvalue of h^hg^- We then choose Wan as the remaining 
— 1 eigenvectors of h^hg,.. Such a design ensures that W 
is a unitary matrix. We then design t as 


t = 


ts 

tAN ’ 


(3) 


where tg denotes the information signal and Ian is an 
{N — 1) X 1 vector of the AN signal. We dehne /3, 0 < /3 < 1, 
as the fraction of the power allocated to the information 
signal. As such, we obtain E [|fsP] = P and E [tANt^t] = 
where E [■] is expectation and Im is the m x m 
identity matrix. Moreover, we confirm that E [xsx|^] = Iat. 
Based on (1), (2), and (3), the received signal at the relay in 
the first time slot is expressed as 


Vr = \J PsdsPh.sr^gtg -f rir, (4) 

where Pg denotes the transmit power at the source, and rir 
denotes the thermal noise at the relay, which is assumed to be 
a zero mean complex Gaussian random variable with variance 
(T^, i.e., rir ~ CN (O, CT^). 

We next express the received signal at a typical eavesdrop¬ 
per located at i, i G $, in the hrst time slot as 

= \JPsd~^hsiWsts + \J Pgd“''hgjWANtAN + n-ii, (5) 

where hg^ denotes the 1 x channel vector from the source 
to the typical eavesdropper located at i, dsi denotes the 
distance between the source and the typical eavesdropper 
located at i, and nn denotes the thermal noise at the typical 
eavesdropper located at i, which is assumed to be a zero 
mean complex Gaussian random variable with variance 
i.e., nn ~ CAf 

In the second time slot, the relay adopts the DF protocol 
to forward signals to the destination. Specihcally, the relay 













first decodes the transmitted signals from the source and then 
broadcasts them after re-encoding. Therefore, we express the 
received signal at the destination as 

Ud = \j Prd~^hrdXR + Ud, (6) 

where Pr denotes the transmit power at the relay, hrd denotes 
the channel coefficient from the relay to the destination, xr 
denotes the transmitted signal of the relay with E [||xr|P] = 1, 
and Ud denotes the thermal noise at the destination, which is 
assumed to be a zero mean complex random variable with 
variance i.e., Ud ~ CAf (O, ct^). 

We next express the received signal in the second time slot 
at the typical eavesdropper located at i is expressed as 


rate of the wiretap code. Re denotes the cost of preventing 
the transmitted wiretap code from eavesdropping, and Pto 
denotes the transmission outage probability (Pto is defined as 
the probability that the instantaneous SNR at the destination 
is less than Tb = 2^’’ — 1). Mathematically, Pto is formulated 
as 

Pto = PT{rD<n). ( 12 ) 

Henceforth, we refer to Tg as the secrecy throughput. 

The goal of this work is to maximize the secrecy throughput 
of the relay wiretap channel with spatially random eaves¬ 
droppers under secrecy constraints. To achieve this goal, we 
formulate the design problem as 


Vi^'’ = \/ Prd^ihriXR + ni 2 , (7) 

where hri denotes the channel coefficient from the relay to 
the typical eavesdropper located at i, dri denotes the distance 
between the relay and the typical eavesdropper located at i, and 
ni 2 denotes the thermal noise at the the typical eavesdropper 
located at i, which is assumed to be a zero mean complex 
random variable with variance cr^ 2 > u,j 2 ~ CAf (O, <^a)- 

B. Received SNRs 

We express the instantaneous SNR at the destination as 


r_D = min{7sr,7rd} , (8) 

where 7 ^^ = and 7 ^^ = For 

r rd d 

the instantaneous SNR at the eavesdroppers, we assume 
that the eavesdroppers are non-colluding, which means that 
each eavesdropper decodes her own received signals from 
the source and the relay, without cooperating with other 
eavesdroppers. Thus, we express the instantaneous SNR at the 
eavesdroppers as' 


where 


Fb = max{max{7sj,7ri}}: 
ie<i> 


PPsdg'’\\hs,ws\ 


and 7 ri = -jr^WKiW" 


N 

C. Problem Formulation 


(9) 

( 10 ) 


In order to evaluate and optimize the secrecy performance 
achieved by our proposed scheme, we apply the performance 
metric proposed in [9], which is given by 

Tg = ^{Rb-Re)(l-Pto), (11) 

where the factor 1/2 is due to the two time slots used in 
the transmission, Re) denotes a parameter pair of the 
wiretap code used by the source, Rt, denotes the transmission 


*We assume that the source and the relay use different wiretap codes with 
different codebooks. As such, the transmitted signals from the source and relay 
cannot be jointly processed by any eavesdropper (nor by any combination of 
eavesdroppers). 


max Tg , 

Rb,Re,fl 

S.t. Pgo < < Re < Rb,fi < P < 1, (13) 

where Pgo denotes the secrecy outage probability. Pgo is 
defined as the probability that F^; is larger than Te = 2 ^® — 1 . 
Mathematically, Pgo is formulated as 

Pgo = Vv{TE>Te). (14) 

III. Analysis and Optimization of Secrecy 
Performance 

In this section, we first analyze the secrecy performance by 
deriving explicit expressions for the transmission probability 
and the secrecy outage probability, respectively. Based on 
these results we determine the optimal parameters, e.g., Rb, 
Re, and /3, that achieve the optimal secrecy performance of the 
relay wiretap channel. Notably, the determined optimal secrecy 
performance is independent of realizations of the main channel 
and the eavesdroppers’ channels. 


A. Transmission Outage Probability 

In this subsection we focus on the transmission outage 
probability, Pto- We first obtain the cumulative distribution 
functions (CDFs) of 7 ^^ and 7 ^^ as 


Fi.r ( 7 ) = 1 - exp 



(15) 


with 7 ^^ = Pgdgf?(7^ and 


Pira ( 7 ) = 1 - exp (16) 

\ 'Ird J 

with 7 ^^ = Prdp^a'^^. Based on ( 8 ), (15), and (16), we re¬ 
express the transmission outage probability in ( 12 ) as 


Pto =Pr (min {7 sr j 7rd} < n) 

=1 - (1 - in)) (1 - (rt,)) 

( 17 ) 



B. Secrecy Outage Probability 

We now focus on the secrecy outage probability, Pgo- To 
derive Pgo, we first express the CDFs of ^si and as 


( 7 ) = 1 - 



(1 -/ 3)7 

/3(iV-l) 


-(AT-l) 


exp 



with = Psdja^^^, and 

F^r-i ( 7 ) = 1 - exp (19) 

with 7 ^^ = Prd~da~ 2 - Based on (9), (18) and (19), we present 
the secrecy outage probability in the following theorem. 

Theorem 1: The secrecy outage of the relay wiretap chan¬ 
nel is derived as 


Fso = 1 - exp (-2A [ Ji + J 2 - J 3 )), 


where 


TT 

PPs 

Jl = — 


77 

\Te(Tf 


J2 = 


-(N-l) 


r( - 


coo pTr 


dsi exp {--ip {6)) ddsidO, 


/o Jo 


( 20 ) 


( 21 ) 

( 22 ) 


Js 



(1 - P ) Te 

PiN-l) 


-{N-l) 


X 



dsi exp 


'TeCrfi 

/3Ps 



exp (—■0 {9 j) ddsidO, 
(23) 


and 0 (0) = {d'^r + dl^ - 2dsrdsi cos 6) ^. 

Proof: See Appendix A. ■ 

We find that (20) provides an easy-to-compute expression 
for the secrecy outage probability. Despite that J 2 and for 
general ry can not be obtained in closed-form, they can be 
easily evaluated since only a double integral is involved in J 2 
and J 3 . 

We next present simplified closed-form expressions for the 

special case where ij = 2. For this special case, we first 

simplify J 2 , which yields a closed-form expression given by 

(24) (next page), where (a) follows by applying the expression 

of the zero-order modified Bessel function Io{z), (b) follows 

by applying the series representation of IqIz), (c) follows by 
2 

using t = dl^ and noting that 

exp (—f) dt = T (k + 1 ) , (25) 

and (d) follows by applying the series expression representa¬ 
tion of the exponential function. Similarly, for the special case 
where ry = 2 , we simplify jTs, which yields a closed-form ex¬ 
pression given by (26) (next page). The simplified expressions 
in (24) and (26) offer us a computationally efficient way to 
calculate the secrecy outage probability for the special case 
where p = 2 . 



C. Throughput Optimization 

In this subsection we determine the optimal parameters 
that maximize the secrecy throughput, Tg, for general ry. 
Specifically, we first determine the optimal wiretap code rates 
pair, {Rl,Rl), that maximizes the secrecy throughput for 
a given power allocation factor (3. We then determine the 
joint optimal power allocation factor and wiretap code rates, 
{l3*°, Rl°, Rl°), that maximizes T^. 

1) Optimal wiretap code rates pair for a given (3: The 
optimal wiretap code rates pair, (i?J, Rf), that maximizes Ts 
for a given (3 is determined as 

{Rl,Rl) = argmax Tg, 
s.t. Pgo < <p,0 < Re < Rb- (27) 

Taking the first-order derivative of Pgo with respect to Re, 
we confirm that dPgo/dRe < 0 , which states that Pgo 
monotonically decreases as Re increases. As such, the value 
of R* satisfying (27) is the value of R* that satisfies the 
secrecy outage probability constraint, i.e., Pgo (i?*) = p. We 
then confirm that dTg/dRt is first positive then negative as 
Rb increases. This demonstrates that the value of satisfying 
(27) is unique. Although a closed-form solution for {Rf Rf) is 
mathematically intractable, we are able to determine the values 
of {Rl,Rl) numerically. Accordingly, the maximal secrecy 
throughput based on the values of {R ^, i?*) for a given [3 is 
defined as T*. 

2) Joint optimization of f3, Rb, and Re-' The joint optimal 
power allocation factor and wiretap code rates which maxi¬ 
mizes Tg in (11), {(3*°, Rpf, Re°), is determined as 

= argmax Tg, 

s.t. Pgo<p, 0<Re<Rb, 0 < /3 < 1. (28) 

Using (17) and (20), we are able to solve (28) numerically. 
Specifically, we first determine the value of using 

(27) for each value of /3. This leads to the secrecy throughput 
with (Rf Rl), denoted by T* = \{Rl- Rl) (1 - Pto). We 
then determine the value of (3 that maximizes T*, denoted 
by (3*°. Accordingly, the value of {Rl,Rf) associated with 
[3*° is determined as {Rl°,Rl°). Finally, the maximal secrecy 
throughput based on the values of /3*°, i?J°, and i?*° is defined 

as t;°. 

IV. Numerical Results 

In this section we present numerical results to validate our 
analysis of the outage probabilities and examine the benefits of 
the proposed scheme. For illustrative purpose, throughout this 
section we concentrate on the practical example of a highly 
shadowed urban area with p = 4. In addition, we adopt Xd^^ = 

Xdl, = 1. 

We first verify the accuracy of the transmission outage 
probability and the secrecy outage probability using Monte 
Carlo simulations. In Fig. 2, we plot Pto versus Tb for 
different values of 7 ^ with = 4 and (3 = 0.5. In this 
figure, we consider 7 ^^ = 7 ^^ = 7 ^. We first see that the 
analytical curves, generated from (17), precisely match the 
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Fig. 2: Pto versus Tb for different values of 7 ^ with 77 = 4, 
TV = 4, and /3 = 0.5. 


simulation points, which demonstrates the correctness of our 
expression for Pto in (17). Second, we see that Pto increases 
monotonically as Tb increases for a given 7 ^, which implies 
that the transmission outage probability increases when the 
transmission rate of the wiretap code increases. We further 
see that Pto decreases as 75 increases for a given Tf,. This 
reveals that the transmission outage probability reduces when 
the source and the relay use more power to transmit under a 
fixed Th¬ 
in Fig. 3, we plot Pso versus Te for different values of 
7 g with = 4 and /3 = 0.5. In this figure we consider 
= lrd^\l<Cl 2 = le- We see an excellent match 
between the analytical curves generated from ( 20 ) and the sim¬ 
ulation points, demonstrating the correctness of our expression 
for Pso in (20). We then see that Pgo decreases monotonically 



T. (dB) 


Fig. 3; Pso versus Te for different values of 7 g with rj = 4, 
TV = 4, and ^ = 0.5. 


as Te increases for a given 7 g, which shows that the secrecy 
outage probability decreases when the redundancy rate of the 
wiretap code increases. We further observe that Pso increases 
as 7 g increases. This is due to the fact the eavesdroppers 
receive signals from both the source and the relay. It follows 
that increasing the transmit power at the source and the relay 
leads to an improved received SNR at the eavesdroppers. 

We now examine the impact of the system parameters Rb 
and P on the secrecy throughput. In Fig. 4, we plot Tg versus 
Rb for different values of N with a fixed /3 and the optimal Re- 
We first observe that there exists a unique that maximizes 
Ts for a given /3. We also observe that the maximal Tg for 
a given /3, i.e., T*, increases as N increases. This shows 
that adding extra transmit antennas at the source signihcantly 
enhances the secrecy performance of the relay wiretap channel. 
































Fig. 4; Ts versus Rb for different values of N with 77 = 4, 
P = 0.5, ip = 0.4, 7 t, = 20 dB, and 7 ^/ 7 ^ = 20. 


relay, we determined the optimal parameters (/3*°, i?*°) 

that maximizes the secrecy throughput, based on our derived 
expressions for the transmission outage probability and the 
secrecy outage probability. In addition, we examined the 
impact of N on the secrecy throughput, showing how the 
maximal secrecy throughput increases with N. The work 
reported here provides valuable insights into the design of new 
physical layer security schemes in which the locations of the 
eavesdroppers are randomly distributed and not known at the 
source. 
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Fig. 5: T* versus /3 for different values of N with p = 4, 
ip = 0.4, 7 f, = 20 dB, and ^bMe = 20. 


Appendix A 
Prooe of Theorem 1 

According to (9), (18) and (19), we re-express (14) as 

Pso = 1 - PrlFfi < Te} 


= 1 - Pr I 

= 1 — Ecj. 

= 1 — Eij. 


max{max{7si,7H}} < 

I 


Pr {max ( 7 ^^, 7 ^} < Te} 


X exp 


Te 


P{N-l) 


hs 


1 — exp —^ 


(a) 


— 1 ~ exp (— 2A (j7i + J 2 ~ iTs)) • 


(29) 


In Fig. 5, we plot T* versus /3 for different values of N. For 
each point of T*, we choose the (PJ, P*) pair that maximizes 
Tg for the corresponding /3. We first observe that there exists a 
unique /3*° that maximizes T*. We then observe that the value 
of P*° that maximizes T* is around 0.9, which demonstrates 
that the secrecy throughput is improved if approximately 10 % 
of the total transmit power at the source is allocated to AN 
signals. We also observe that the maximal T*, i.e., T*°, 
increases as N increases. Furthermore, we observe that the 
value of P*° slightly decreases as N increases, which shows 
that in order to maintain the optimal secrecy throughput, the 
power allocated to AN must increase as the source antenna 
number increases. 


V. Conclusion 

In this work we proposed a secure transmission scheme for a 
relay wiretap channel, in which an N-antenna source transmits 
both information signals and AN signals in the presence of 
multiple spatially random single-antenna eavesdroppers. Con¬ 
ditioned on the use of a decode-and-forward protocol at the 
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r„cr. 


PPs 


dl-^dP]dds^d0, (32) 


and (a) follows by applying the probability generating func¬ 
tional (PGFL) for the PPP $, given by [30] 


E$ 


n/i 


= exp{-yjl-/(a;)]Afia;}, (33) 


















and by changing to polar coordinates. 
To proceed, we first derive Ji as 
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(34) 


where in (b) we have used u = in (c) we have used 
t = ^ and (d) follows from the definition of the gamma 

function. We then derive Ji as 


Ji 


dsi exp - 


0 JO 


Pr- 


'-{dl^ + d4-2dsrdsi cosd) ^ ^ddsidO. 

(35) 


We further derive as 


Jo 



Substituting (34), (35), and (36) into (29), we obtain the 
desired result in (20), which completes the proof. 
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